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Abstract 

Some magnetoelectric multiferroics have a canted spin structure that can be described by a 
Dzyaloshinkii-Moriya coupling. We calculate properties and features expected for surface and bulk 
magnon polaritons in such media with a linear magnetoelectric interaction for the case of transverse 
electric polarisation. The dielectric polarisation and magnetisation of weak ferromagnetism are 
constrained to lie in the plane parallel to the surface. We examine a geometry with the polarisation 
oriented in the film plane and present numerical results for the transverse electric polarisation. 
Particular attention is given to non-reciprocal surface modes, which exist in frequency between 
two bulk bands, and show how these modes can be modified by external magnetic field. Results 
for attenuated total reflection are presented, and discussed in relation to nonreciprocity. Example 
results are calculated for the canted antiferromagnet BaMnF4. 

PACS numbers: 71.36+c; 78.20.Jq; 78.20.Ls 
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I. INTRODUCTION 



Polaritons in magnetic and dielectric media arise from coupling between the electromag- 
netic field and polarisations of the mediajlj]. A particularly interesting polariton is localised 
to the surface of a bounded material. Surface polaritons can display nonreciprocity, where 
the propagation frequency in one direction is different from the propagation in the opposite 
direction. If the propagation wavevector is k, then nonreciprocity can be expressed in terms 
of frequency wasw ^ to [— k). Numerous applications for nonreciprocal surface waves 



exist [2|, |3[ including recent developments for surface polariton optics 



Magnon po 



ferromagnets 



aritons are coupled photon and magnetic excitations, and_exist for simple 

In the present 



, |7| and multisublattice magnets including antiferromagnets 
work, we discuss surface modes for an interesting class of materials that are currently a focus 
of attention: magnetic multiferroics. Multiferroics display long range order and polarisable 
response in two or more aspects: elastic, dielectric and magnetic. The calculations pre- 
sented here are illustrated using parameters appropriate for BaMnF4, a material that has 
been modeled with a linear magnetoelectric coupling between the magnetic and dielectric 
subsystems. This material has two magnetic sublattices, and a canting angle between sub- 
lattice magnetisations results from magnetoelectric interaction {9 -11]. 

Some properties of bulk polaritons in linear magnetoelectric coupled media have been 
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and surface modes have been examined in Refs. 



culated theoretically in Refs. |9| and 

and Ji! To the best of our knowledge, until now dispersion relations in the presence of 
applied electric and magnetic fields have not been discussed, and explicit results for surface 
modes in BaMnF 4 have not been presented. 

In this work, we show that nonreciprocity of magnetic polariton surface modes can be 
modified by application of an external magnetic field. We extend previous work by allow- 
ing magnetic sub-lattices to cant, and provide explicit results for surface and bulk modes 
on BaMnF4. Our results are relevant for understanding and predicting microwave and in- 
frared responses of ferroelectric/magnetic and multiferroic/magnetic multilayers fl^. [lol] and 
composites Q, Q- In these types of heterostructures, it may be possible to achieve effective 
media with magnetoelectric properties and a range of material parameter values not possible 
in single phase multiferroics. 

The paper is organised as follows. The field dependence of the magnetic canting angle 
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is discussed in the next section, and relevant magnetic and dielectric susceptibilities are 
calculated in section III. Results for surface and bulk mode dispersions are presented and 
discussed in the final section IV. 

II. GEOMETRY AND MAGNETIC SUB-LATTICE CANTING 




FIG. 1. Geometry. Canting of two magnetic sublattices (m a and nib) by an angle 9 produces a 
weak ferromagnetism (M) along the x axis parallel to the surface. The spontaneous polarization 
(P) is assumed to lie in planes parallel to the surface and propagation of the surface mode is along 
the y axis with wavenumber ky. 

The geometry used is shown in Fig. [TJ A semi-infinite film of a multiferroic with a two 
sub-lattice antiferromagnet fills the half space z < while the other half space is assumed 
to be vacuum. The magnetocrystalline easy axis for a uniaxial anisotropy is directed out of 
plane along the z axis. The spin system cants symmetrically in the x — z plane, and the 
sublattices are assumed to have the same magnitude of magnetisation, \m a \ = |m&| = M s . 
The nett magnetisation, M corresponding to the weak ferromagnet due to spin canting, lies 
along the x axis in a plane parallel to the surface. The spontaneous dielectric polarisation is 
oriented in the perpendicular direction, also in the x — z plane. An external magnetic field 
is applied parallel to the weak ferromagnetic moment. 

The transverse electric (TE) mode is associated with the condition where the electric part 
of electromagnetic waves propagate along the surface is parallel to the surface. We assume 
that the surface modes propagate on the surface along the y direction, hence the electric 
components lie on the surface in the x direction while the magnetic part has H y and H z 
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components. 

A fourth order Landau- Ginzburg density energy is assumed for the spontaneous ferro- 
electric polarisation: 

Fe = ^ClP* + \<tf* - PyE y . (1) 

The coefficient Ci is temperature dependent and an external static applied electric field E y 
along spontaneous polarisation is included. 

The magnetic contribution to the density energy is given by: 

K 

F M = -Xrh a ■ rh b - — - [{rh a ■ zf + {rh b ■ z) 2 ] - {rh a + rh b ) x H (2) 

where A and K represent magnetic exchange and anisotropy energies, rhi is the magnetisation 
of sub-lattice i and H Q is an external static magnetic field. The first and second terms in the 
right hand side Eq. (j2J) represent the exchange and anisotropy energy, while the third term 
is Zeeman energy contributed by external magnetic field. In this configuration, the easy axis 
is out of plane. We assume that there is zero nett magnetisation along the z direction and 
do not consider demagnetisation fields. 

The magneto-electric coupling is assumed to be of the Dzyaloshinkii-Moriya (DM) form, 
appropriate to materials such as FeTiC>3 18] and BaMnF4js-.ll]. We suppose a linear cou- 
pling of the form 

F me— -aP y M x L z 



-aP y ■ (rh a x fh h ) - aP y [{m a ) x (m a ) z - (m b ) x (m 6 )J (3) 



where 



M x = {m a + rh b ) x = 2M S sin 9, (4) 

L z = {rh a - rh b ) z = 2M S cosd (5) 

and a is the magneto-electric coupling constant. The component M x is the transverse 
component of magnetisation and represents the weak ferromagnetism. The component L z 
is the magnitude of the longitudinal component of the magnetisation. The first term in Eq. 
(J3]) represents the polarisation induced Dzyaloshinkii-Moriya interaction in our model, and 
we assume that this term is responsible for generating canting in the magnetic sublattices. 

The canting angle is calculated by minimizing the magnetic and magnetoelectric free 
energies. Writing {F M + F ME ) = 0, we obtain the condition: 

H cos 9 - ^KM S sin 26 + 2aP y M s cos 29 + \M S sin 29 = 0. (6) 



This can be written in the form 

4aP 

tan (7) 

when the external magnetic field is absent, H a = 0. As can be seen from Eq. ((Zj) above, the 
magnitude of the magnetoelectric coupling and spontaneous polarisation control the canting 
angle. 

The equilbrium spontaneous polarisation is found by minimizing the dielectric and mag- 
netoelectric free energies, -^p (F E + F ME ) = 0, with respect to the spontaneous polarisation, 
with the result: 

&P V + &Py - 2«M S 2 sin 26 - E y = 0. (8) 

Lastly, the temperature dependent magnetisation in mean field is calculated using the Bril- 
louin function in B (if) 

M s (T)=M s (0)B s ( V ) (9) 



where 



= 9^bS_ = r_ 2Q + R cog2 e + 2aP M gin 29 + Ho gin q\ _ / 10 n 

k B T L J 



The spontaneous polarization, magnetization and canting angle are obtained by solving Eqs. 
(EJ), (|H]) and (jHl) simultaneously using numerical root finding techniques. 

III. DYNAMIC SUSCEPTIBILITIES 

The electromagnetic problem requires constitutive relations for magnetisations and po- 
larisations. These are formulated in terms of linear magnetic and dielectric suceptibilities 
and are derived from the magnetic Bloch equations: 

M = 7 x (^(Fm + Fme)) (11) 



and the Landau-Khalatnikov dynamic equations for the polarisation P: 

P = -fJL(F E + F ME ). (12) 
oP 

Here 7 is the gyromagnetic ratio and / is the inverse phonon effective mass. 
The relevant equations of motion for TE modes after linearisation are, 

— iul x = 2 / ~fM s hy cos 9 — (u a cos 6 + 2u me sin 6 — 2u ex cos 9) m y , (13) 
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— iurniy = 2 r yM s h z sin 9 + (u a cos 9 + 2oj me sin 9) l x + (u a sin — 2u me sin 6 1 — w a ) m 2 , (14) 
— io;m 2 = 2 r yM s h y sin 6 1 + (u; + 2u me cos 6 1 ) m y (15) 

and 

^y-P* = ~ (Ci + C2P0) p a + e x (16) 
where Zj = m ati — m b i and 77^ = m a ^ + 77^. The frequencies u a = jKM s , u ex = 7AM,, 
w me = 7aP G M s and co = •yHo are associated with the magnetic anisotropy, exchange, 
magnetoelectric coupling and the external magnetic field respectively. The parameter P Q 
represents the equilibrium spontaneous polarisation. It can be seen from Eqs.f lT3|) to ffl6l) . 
that the dynamic magnetic and electric components are not coupled directly. Hence, the 
susceptibilities can be written as 

rh = x m h and p = x e e, (17) 

with magnetic components given by 

m _ 2jM s (uj a cos 29 + 2u me sin 29 + u sin 9) 

Xy ~ H fm cos*9 + ni e + ni-ui) [ } 

m _ 2jM s (ay e sin 29 + u D sin 9) 

Xz ~ H fm ^9 + ni e + ni-^) [ } 

XyZ " XzV K m cos^ + fiL + ^-^ 2 ) 1 } 

Here = co a (cj a — 2u ex ) represents the antiferromagnet (AFM) resonance frequency. The 

frequency Q me is defined as 

,2 - /- 1 



VL me = 2u me \ 2u me + -u a sin 29 - u ex sin 29 j (21) 

The frequency Q Q is the frequency shift due to an external magnetic field and is given by 

0% = co [tu — co a sin 9 + Auj me cos 9] (22) 

Even though there is no direct coupling between magnetic and electric parts represented 
by a ME susceptibility, the ME u me still appears in the magnetic susceptibility, as seen in 
Eq. ( fl8l) to ( |20l) . This results from the fact that the canting condition involves the electric 
polarisation. Finally, the electric susceptibility for \ e xx 1S 

X-tx = 77§ ( 23 ) 



where ujl is the frequency of the longitudinal phonon mode, and uy is the frequency of the 
transverse phonon mode. 



IV. RESULTS AND DISCUSSION 

Since there is no dynamic ME coupling, the dispersion relation for bulk polaritons can 
be calculated from the macroscopic electromagnetic wave equation 

V 2 H - V (V • itj - ■ d 2 H/dt 2 = 0. (24) 

For the TE modes, the constitutive connecting B with H are defined by jl = 1 + 47rx m . 
Using these relations, the solution of the wave equations for the bulk modes is of the form 

H ~ e i(k y y+k z z-ujt) ^25) 

and an implicit expression for the bulk mode frequency is: 

Ihtfi = e * (~) (/V* ~ AV 2 ) • ( 26 ) 

Equation (1261) has two zeros, one is for e x = and the other for / (/i) = fi y fi z — fiy Z = 0. It 
also has two resonance poles from e x and another from the condition fi y = 0. As a result, 
there are three bands of bulk polariton modes. 

The dispersion relation for surface modes is obtained by assuming solutions in the form: 

H ~ e ^ e *(M-^) for z < (27) 

and 

H ~ e -Poz e i(k y y-ut) f or 2 > (28) 

where /3 and (3 are positive real attenuation constants for the sample and vacuum, re- 
spectively. An implicit relation for the attenuation factor (3 of the medium is derived by 
substituting Eq. (j27|) into the wave equation Eq. ff24l : 

fi z P 2 = fi y k 2 y - e x (^j (fi y fi z - fi 2 yz ) (29) 

An explicit relation for the attenuation constant /3 Q (in vacuum) is given by 
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K = K-{-)- (3°) 
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An implicit surface mode dispersion relation is calculated by requiring continuity of tangen- 
tial H and normal B at the interface z — 0. Using Eq. fj29|) and ( 130]) . we find 



Hzfl + {j^y^z l^yz) 0O "F Hyzky 



0. 



(31) 



The dispersion implied by Eq. (l3Tj) involves mainly of magnetic permeabilities, and so 
describes magnetic surface polaritons for the weak ferromagnet. From symmetry, we know 
that magnetic surface polaritons can be nonreciprocal 19(, and this is the case here as well. 
The direction of k y matters in Eq. ( 13T1) and indicates nonreciprocal propagation. 

A key point is that the nonreciprocity of the surface modes depends strongly on the 
canting angle. If the canting angle is zero, then fi yz — > and ji z — > 1, and the dispersion 
relation in Eq. ( )3T|) becomes 

+ VyPo = (32) 

illustrating the dispersion of the surface modes in this case is reciprocal under reversal of k, 
and propagation of the surface modes is reciprocal. 
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CO/27Tc (cm l ) 

(a)Attenuation constant. 
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k (cm -1 ) 

(b)Dispersion relation. 



FIG. 2. Attenuation constant and Dispersion relation. In (a) the attenuation constant is shown for 
two values of wavevector in the absence of external fields. The solid line represents k = cm , 
while the dashed line correspond to k = 16 cm" 1 . In (b) bulk and surface mode dispersions are 
shown. Surface modes are indicated by SP. The shaded regions represent bulk bands, limited by 
frequencies cv m , lo p and w 2 . The dashed lines denote by vertical and horizontal arrows represent 
ATR light lines for incident angles of 30° and 70°. The asymptotic boundaries and lightline are 
represented by dashed lines which are indicated as A and by L. 
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Solutions of Eq. fl26|) and Eq. ( I31j) are plotted in Fig. 2(b) for the case of no applied fields. 
Before continuing, we discuss our choice of parameters, which were chosen to be appropriate 
for BaMnF 4 . The sublattice magnetisation is approximated by M r = 2M S sin 9. Assuming 
the weak ferromagnetisation, as mention in review paper by Scott[20j], M c = 1460 A/m with 
a canting angle of 3 mrad at 4.2 K as reported by Venturini 2l|, results in M s = 3.054 x 10 3 
Oe. The exchange field is estimated by Holmes 22|, as He=50 T. Using the definition He = 
—XM S , we obtain A= -163.72. Since antiferromagnet resonance was reported 23| at around 
3cm -1 , the anisotropy constant is estimated by using expression w 2 = j 2 K(K — 2A)M 2 , 
with the result that K = 0.337. Here, we estimate the gyromagnetic ratio by calculating 
7 — Q^s/h, yielding 7 = 0.933 cm _1 T -1 . 

The spontaneous ME coupling a is calculated from Eq. (j7j) using the calculated sponta- 



neous polarisation 



241 ]. P o =0.115 C/m 2 , from which we obtain a— 4.27 m 2 /C. The phonon 



inverse mass is calculated using the relation between phonon frequency in transversal and 
longitudinal modes, u L = u^ + 47r,f, yielding / =6.4458 cm~ 2 (where ojl = 41 cm -1 and 
ujt = 40 cm -1 as given in Ref. [23J and 20J]). The Ginzburg-Landau constants (j and 
(2 are calculated by solving Eq. (JHJ) with the equation for the soft phonon frequency, 
(Ci + X2P0) f = wii resulting in & = -1.22 x 10 2 and C2 = 1-085 x 10~ 7 cm 4 /statC 2 . 

We also note that previous investigations 23J] have assumed the dielectric constant in the 
x direction to be independent of frequency with the value e x = 8.3. This assumption is valid 
if the dielectric and magnetic responses lie in very different frequency ranges. As a result, 
only two bands for bulk polaritons can exist since only one pole and one zero are contributed 
by the magnetic sub-system. We make this assumption in what follows, since BaMnF4 has 
a wide separation in frequency between the dielectric and magnetic resonances. At the end 
of this section we discuss the consequences of relaxing this assumption. 

The bulk mode frequencies are affected by the attenuation constant, since this parameter 
determines the frequency and wavenumber regions that allow surface modes. A positive 
value of real (3 is needed for the solution of surface modes in Eq. ( 1271) . Using the above 
parameters, the bulk mode region is bounded by k = cm" 1 and k=16 cm" 1 as shown in 



Fig. 2(a), and is indicated in Fig. 2(b) by shading. Comparison of Figs. 2(a) and 2(b) 



shows that the resonance at 0J m = 3 cm" 1 is associated with a divergence of the attenuation 
constant due to the zero of \i z . From this zero condition, an expression for u m can be 
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derived : 

= [(co 2 afm cos 2 9 + n 2 me ) + (2u s u me sin 29)] 1/2 (33) 
with u s = ^AttM s . 

The pole and zero frequency of bulk modes, u p and u z , is related to the zeroes of fi y and 
f(n) = Hyfi z - fx 2 yz . These give 



u„ 



[(u 2 afm cos 2 9 + Q 2 me ) + {2u s u a cos 29 



1/2 



and 



where 



U, 



1 

71 



tt 2 mp + (tt 4 mp -4u 2 m u 2 ) 



2\V2 



1/2 



Q 2 

mp 



u J 2 m + u 2 p + Au; 2 sm 2 9 



(34) 



(35) 



(36) 



These three frequencies, u m , u p and u z divide the bulk region into the shaded regions shown 



in Fig 2(b) Note that the lowest two bands overlap, and are deformed strongly near the 
resonance u m . 



Surface modes, which is indicated by "SP" in Fig. 2(b) can exist in the gaps between 
bj m and Up. The surface branches are nonreciprocal with respect to propagation, with 
u (jtj 7^ u(—k). The negative branch begins from the intersection of the lightline with the 
u m resonance, and is indicated by "-L" . The branch terminates in the left middle bulk 
band. The positive branch begins at the lightline with frequency 

~£xui-u 2 afm cos 2 9 + n 2 m 1 1 - 



u 



(37) 



e x - 1 

which is slightly higher than u m , and then terminates at the upper bulk band. 

A method for detecting surface modes at THz frequencies is Attenuated Total Reflection 
(ATR)[l] and has been used to study surface modes on antiferromagnets 25]. In this method, 
a high optical index prism is used to couple electromagnetic radiation to surface excitations 
that lie off the vacuum light line. A sharp dip of reflectivity indicates excitation of surface 
modes. Using results from Ref. js], reflectivity in ATR is given by 

k z (1 + re- 2l3 ° d ) - i/3 (l - re' 2 ^) 



R 



k z (1 + re- 2 P° d ) + i/3 Q (1 - re~ 2 ^ d ) 

where d represents the distance between prism and the sample, and r is defined as 

{3 - K 
r = 

Po + K 



(38) 



(39) 
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with 



Here, the wave- vector 



and 



HyHz ~ p 



IF- 



k v = (e n ) 1/2 -sm6 



h = (en) 1 ' 2 -cos9 



(40) 
(41) 
(42) 



represents propagation along and normal to the surface, where e p is the prism dielectric 
constant, and 9 is the incident angle. 



^0.5 




CQ/27Cc(cm L ) 

(a)ATR spectra at 30° 



CO/27lc(cm j 

(b)ATR spectra at 70°. 



FIG. 3. ATR spectra with incident angles 30° and 70°. In (a) two different sharp dips illustrate 
the non-reciprocity of the surface modes. In (b) the absence of the sharp dips indicate the absence 
of surface modes. 



Calculated ATR results for the BaMnF 4 is presented in Fig. 3(a) and 3(b) for the incident 



angles 30° and 70°. The ATR spectra at 30° illustrates nicely the nonreciprocity of surface 
modes. Two sharp dips correspond to surface polaritons traveling in the opposite directions 
and the frequencies are the u m intersection points discussed above for the positive and 



negative surface mode branches. The ATR for an incident angle of 70°, shown in Fig 3(b) 
does not allow coupling to surface modes. The shallow dips represent bulk modes. 

It is interesting to notice that the surface modes do not exist in the region where the 
wave- vector k 3> -. In this region, the surface dispersion relation in Eq. ( 13 ip can be 
re-written in the form 



fiyfiz ~ {Hyz =F I) 2 = 



(43) 
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for k = ±00 . In terms of u m , u p and u r Eq. (|43|) becomes 



( 



)(<■ 



UJ m — US I I L0„ — LO 



1/2 

where uj r = {oj 2 a ^ m cos 2 6 + fi^e) represents the pole frequency of permeabilities given 
by Eq. (ll3H16]) and is located slightly below uj m so that u r < u m < u p . Assume that the 
asymptotic frequency for the surface is located in the region above the pole frequency u p . 
For k = +00, the second term in Eq. (jUJ) is always larger than the first term. Hence, Eq. 
is never satisfied and the surface modes do not exist in this wavelength region. 



2 ) — [2u; s u; sin 9 =1= u 



.2M2 



)] =0 



(44) 
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(a)The absence of surface modes at —k. (b)The existence of surface modes at —k. 
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(c)The parameters are A/a=1.5 



k (cm ) 



(d)The parameters are X/a=2. 



FIG. 4. The function in Eq. (|43p represented by f(n) is shown in (a) using parameters for BaMnF4. 
The condition for surface modes in k >> u/c region has no solution. In (b) a solution exists using 
parameters where A/a = 1.5, where the surface begin to overcome the bulk as shown in (c). In (d) 
the surface branch is slightly higher than bulk when the ratio increase to A/a = 2. 



A surface mode at k — —00 may exist for some values of u r , u m and u p . For a solution to 
exist in this region, the function /(//) of Eq. f )44|) should cross zero. This gives a condition 

12 



from Eq. ( 1441 at the frequency u> = u p for a surface mode solution to exist. The requirement 
that / (/i) vanish means that 



u p < uj a cos 29 + 2u me sin 29 



(45) 



For BaMnF4, the requirement in Eq. fj43j) is never satisfied as shown in Fig. 4(a) , hence 



there is no surface modes at k » —. The condition Eq. ( I4"5j) can be satisfied by changing, 



for example, the parameters so that A/a = 1.5. An example is shown in Fig. 4(b) , and the 



associated dispersion relation is presented in Fig, 4(c) In this case a surface mode is able to 



exist at the top of the middle bulk band, for arbitrarily small wavelengths. 
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(a)canting angle versus H 
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(b) Surface modes at H = 15T 



FIG. 5. The influence of magnetic field. In (a)influence on canting angle. In (b), the application 
of a magnetic external field of 15 T increases and allows the negative surface branch to enter the 
magnetostatic region, where k >> — . 



In principle, application of an external field can be used to modify the surface mode 
frequencies by changing the canting angle. However, as mention above that the electric part 
is uncoupled from the magnetic part, and application of an electric field has negligible effect 
on the polariton frequencies. 

The application of an external static magnetic field changes directly the canting angle 
and thereby affect the frequencies. The canting angle as a function of applied magnetic field 



shown in Fig. 5(a) In fact, for sufficiently strong magnetic field, the condition of Eq. fj43|) 
can be satisfied and a surface modes branch extend to k = — oo. This occurs through an 
increase in u r via the field dependence of n o (see Eq. (|22l)). 
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The frequencies u m and pole frequency u p also increase under an external magnetic field. 
The difference between uj m and u p is given approximately by 



A 



uj 2 m ) oc (cos 29 - sin 29) 



(46) 



As a consequence, the middle bulk band narrows in frequency with increasing magnetic field. 
At a field of 15 T, the BaMnF4 middle bulk band is sufficiently narrow that a surface mode 



rises above the middle band and extends to small wavelengths, as illustrated in Fig. 5(b) 

Lastly, we consider the case where the dielectric constant depends on frequency. We 
suppose that the frequency dependence is given by the soft phonon mode according to 



c4= ((3 1 + P*f3 2 )f 



with the dielectric function 



(47) 



(48) 



Here u;£ = ut. + — and e& = 8.3 is dielectric constant background. 
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(a)The case where w e > ui m . (b)The case where uj rn > uj e . 

FIG. 6. Dispersion relation for a frequency dependent dielectric function. The frequency is in- 
troduced according to Eq. (35). In (a), the case when of a magnetic resonance frequency lower 
than the electric resonance frequency is shown. In (b), the case of a magnetic resonance frequency 
higher than the electric resonance frequency is shown. 



Since there is no dynamic electric and magnetic coupling for the TE modes, there is no 
effect on the surface mode frequencies even though the electric resonance is near the magnetic 
resonance in frequency (this can be achieved by modifying the phonon mass). Effects on 
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the bulk bands for an inverse of the phonon mass of 0.015 / is shown in Fig. 6(a) where 
the magnetic resonance is lower than the dielectric resonance. The effects of a magnetic 



resonance above the dielectric resonance, with 0.009/, are shown in Fig. 6(b) Only band 
overlap and intersections are affected. Note a consequence of this for the surface modes. 
In the case where the magnetic resonance is below the dielectric resonance, have a broader 



range in wavevector (see Fig. 6(a)) compared to the case where the magnetic resonance is 



above the dielectric resonance (in Fig. 6(b)). 



V. CONCLUSIONS 



Dispersion relations for a material having linear magnetoelectric coupling have been cal- 
culated for tranverse electric polariton modes. Effects associated with spin canting are 
considered. The polarisation and weak ferromagnetism are parallel to the surface, and the 
antiferromagnet magnetisations are out of plane. We find surface modes associated with the 
weak ferromagnetism that are non-reciprocal with respect to propagation direction, such 
that tu (jtj 7^ co(—k). For sufficiently large magnetic fields, or different material parameters, 
surface modes may also exist in the magnetostatic region, with k ^> u/c. Application of 
static magnetic field can be used to modify the middle bulk band frequencies 
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